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ABSTRACT 

The expanding ejecta from a classical nova remains hot enough 10"^ K) to be detected in thermal 
radio emission for up to years after the cessation of mass loss triggered by a thermonuclear instabil¬ 
ity on the underlying white dwarf (WD). Nebular spectroscopy of nova remnants confirms the hot 
temperatures observed in radio observations. During this same period, the unstable thermonuclear 
burning transitions to a prolonged period of stable burning of the remnant hydrogen-rich envelope, 
causing the WD to become, temporarily, a super-soft X-ray source. We show that photoionization 
heating of the expanding ejecta by the hot WD maintains the observed nearly constant temperature of 
(1 — 4) X 10^ K for up to a year before an eventual decline in temperature due to either the cessation of 
the supersoft phase or the onset of a predominantly adiabatic expansion. We simulate the expanding 
ejecta using a one-zone model as well as the Cloudy spectral synthesis code, both incorporating the 
time-dependent WD effective temperatures for a range of masses from 0.60 Mq to 1.10 Mq. We show 
that the duration of the nearly isothermal phase depends most strongly on the velocity and mass of 
the ejecta and that the ejecta temperature depends on the WD’s effective temperature, and hence its 
mass. 

Subject headings: stars: novae - stars: cataclysmic variables - stars: white dwarfs - X-rays: binaries 


1. INTRODUCTION 

Classical novae (CNe) are caused by the ejection of 
AMej ~ 10“^ Mq of matter from an accreting white 
dwarf (WD) triggered by a thermonucle ar instability 
in recently-accreted hydrogen-rich ma tter (Gallagher & 
Starrfield 1978| |Starrfield et al.||2008|). The expanding 
ejecta are bright m the radio tor hundred s of days after 
the optical peak (Seaquist & Bode 2008) and are well 
modeled by assuming the radiation comes from thermal 
bremsstrahlung emission from an expanding medium of 
nearly constant electron temperature Tp » (1 —4) x 10^ K 


dSeaquist fc Palimaka|1977 Hiellming et al.|I979 Taylor 


et al.||1988| |Eyres et al.||1996|). buch temperatures are 

also observed via nebular spectroscopy lo ng after opti¬ 
cal p eak in novae like V1974 Cygni 1992 ( Aust in et al. 
1996|) and T Pyxidis V35I (|Shore et al.||2U13|). How- 
ever, between the fading of the optical einission (which 
may correspond to the end of mass ejection) and the 
end of the radio emission several hundred days later, the 
ejecta undergo an expansion in volume by several orders 
of magnitude. In the absence of a heating source, the 
ejecta would undergo adiabatic expansion, leading to a 
large temperature decline to values orders of magnitude 
less than the observed 10'^ K. The cause of these much 
higher temperatures for such prolonged periods is the 
focus of this paper. 

With the rise of frequent X-ray monitoring of classi¬ 
cal novae, it is becoming increasingly clear that most, if 
not all, CNe are followed by a super-soft X-ray source 
(SSS) phase during which the WD shrinks in radius to 
nearly its pre-nova size, but continues to burn a rem¬ 
nant hydrogen envelope at nearly the Eddington lumi- 
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nositv as modeled by Starrfield et al. (II974), Fujimoto 

(1982), Sala & Hernaii 

z (2005), Hachisu fc Rato (|20iO), 

and Wolt et al. 121)13) 

This phase has been decteced in 

many galactic 

lovae (Uegeiman et al. II 1984 

Shore 2008 

Krautter|2008 

Schwarz et al.||2Ull), and obseryations in 


follow every classical nova outburst (Henze et al. 2011 
2014bI. During the SSS phase, the WD has L « 10 Lq 
and since it is compact, is a strong source of radiation at 
fcTeff « 50 — 100 eV. This SSS phase can be as short as 


a few days for novae with the m ost massive WDs (Tang 
et al.|[2014 Henze et al. 2014a) a nd as long a s a decade 


tor low-mass WDs like V723 Gas (Ness et al. 


2008). We 


show here that the SSS phase is a sufficiently powerful 
source of photoionizing heating to sustain the high ejecta 
temperature long after the optical peak. 

We start in Section by reviewing the physics of pho¬ 
toionization balance and establish our simplified model 
of hydrogen and helium ionization in nova ejecta. We 
derive the temporal evolution of the ejecta temperature 
in Section and discuss computational models in Section 
1^ We compare our semi-analytic results to the computa¬ 
tional results in Section]^ and end by showing how this 
new understanding can inform us about the properties 
of CNe ejecta, and, potentially, the WD mass. 


2. PHOTOIONIZATION BALANCE 

For the late time evolution of interest to us, we chose 
to model the nova ejecta at age t as one zone of constant 
mass AMej in a uniform density sphere (i.e. not a shell) 
of radius r{t) = v^jt expanding at constant velocity Uej. 
A hot WD of mass M and radius i?wD resides at the 
center. The ejecta has a composition X, here represented 
as a vector of mass fractions. The number density of any 
given isotope, i, in the ejecta is 


n(W) 


3WAMej 

AiUlpATTV^-fi ’ 


( 1 ) 
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where Ai is the mass number of isotope i and rUp is the 
proton mass. The electron-ion coupling time is always 
much shorter than the age, t, during the time relevant for 
our calculations, so there is only one temperature, T{t). 
We express the fundamental equation of thermodynamics 
in terms of ejecta volume V, entropy S, and total number 
of particles N as 


TdS = NkT 


[MT 



( 2 ) 


For a spherical ejecta of pure ionized hydrogen {N = 
2Ne) with r = vt, this becomes 


T 


ds 

dt 


UpkT 


S^lnT- 

dt 


(3) 


where we now consider the specific entropy s = S/V and 
rig is the number density of electrons. 

The only source of heating we will consider is that from 


tern that has been modeled extensively dStarrfield et al. 

19741 Fuiimoto|1982l MacDonald & Vennes|1991l 


Hernanz 2005[ IHachisu & Katol 20101 

Wolf et al. 

2013). 

We extended the work ot Wolf et al.j ( 

2U13) to generate 


as shown in Figure These simulations used the 1-D 
stellar evolution code MESA (|Paxton et al. 2011 2013), 
following the time-dependent evolution ot wDs accret¬ 
ing solar composition material at 10“® Mq yr“^ through 
several nova flashes with super-Eddington winds as the 
primary mode of mass loss. In particular, we only dis¬ 
play the results from the end of mass loss until the end 
of quasi-stable hydrogen burning for one representative 
nova cycle. The original simulations did not have data 
for a 0.8 or 1.1 Mq WD, so we generated additional data 
for this work. At any given age since mass loss ended, 
we know the WD radius and Teffi and hence the radia¬ 
tion field that powers the photoionization heating. This 
is the only information we use from the MESA models. 
The ejecta parameters AMej and 'Cej used are indepen¬ 
dent of the MESA models, and, in general, would vary 
significantly for novae on WDs of different masses. 


2.1. Stromgren Breakout 

Our first concern is whether the photoionization source 
is adequate to maintain a large region of mostly ion¬ 
ized plasma in the ejecta. Hence, we start by finding 
the Stromgren sphere radius, Rg, and compare it to the 
ejecta dimension, Vejt. For simplicity, we assume pure 
ionized hydrogen (rie = rip) within the Stromgren sphere, 
and equate the rate of emission of H-ionizing photons by 
the WD, Qn, with the rate of recombinations in that 
volume, 

Qh = -nRlanpUe , (4) 

where a is the radiative recombination rate of hydrogen. 

At early times Rg is far inside the outer ejecta radius, 
Uejt, and the ratio scales as 


V^jt 


/ dTrQn 

\3aN^ 


1/3 


{Vejt), 


(5) 


where Np is the total number of protons. Hence, at early 
times, the Stromgren sphere is inside the ejecta, and only 


Table 1 

Stromgren Breakout Times'^ 




H 


He 

Mwd 


Qh,47‘^ 


QHe,47® 

(Mq) 

(d) 

(10"^'^ s-i) 

(d) 

(10"*^ s-1) 

0.60 

31 

8.0 

95 

1.2 

0.80 

28 

11.0 

69 

3.6 

1.00 

31 

8.0 

60 

5.6 

1.10 

37 

5.0 

67 

4.4 

1.20 

41 

4.0 

214 

3.6 

1.30 

132 

2.9 

387 

2.7 


^ AMej = 10 ® Mq; Uej = 10® kms ® 

Stromgren breakout time for pure hydrogen 
^ Emission rate of H-ionizing photons 
Stromgren breakout time for pure helium 
® Emission rate of He-ionizing photons 

breaks out at an age of 


ts ~ 30 d 


10® km/s 


AMej 
10-5 Mq 


2/3 


Qh 

10^8 S -1 


- 1/3 


(6) 

where we have scaled Uej, AMej, and Qh to values typi¬ 
cal of novae on intermediate-mass WDs. Table [T] shows 
that the H Stromgren sphere breaks out by 50 days for 
all but the most massive WDs for AMej = 10“® Mq 
and Vej = 10^ kms“^. This confirms a time after which 
all of the ejecta is exposed to the photoionizing radia¬ 
tion field. The pure He breakout times for singly ion¬ 
ized He Stromgren spheres are systematically later due 
to less ionizing radiation available for a pure He ejecta. 
The more massive WDs have longer Stromgren break¬ 
out times at the same ejecta parameters because while 
they produce harder spectra than lower-mass models, the 
total luminosity is not significantly higher, so the to¬ 
tal number of ionizing photons emitted, Q, decreases, 
requiring more time for a complete ionization of the 
ejecta. Additionally, the SSS phases for the most mas¬ 
sive models are short enough that the ionization front 
never makes it to the ejecta radius before the ionizing 
source fades greatly, causing excessively long Stromgren 
breakout times. We explore the ionization state of the 
ej ecta at late t imes i n Section!^ 


Beck et al.| (1990) performs a much more detailed 
calculation oi this breakout for a He Stromgren sphere 
for a mixed ejecta with heavier elements {Z > 2), and a 
lower ejected mass and velocity (AMgj = 1.3 x 10“® Mq 
and Uej = 6 X 10^ km s“^ than we used for Table So as 
to make a close comparison, we used their parameters , 
including their WD Teg (see Figure 6 of Beck et al.|1990| 
as inputs to our homogenous model. We found that the 
H and He Stromgren breakouts occur at 31 days and 
77 days, respectively, very close to the breakout they 
reported at f = 76 d for the He Stromgren sphere. 


2.2. Ionization Timescale 

After the Stromgren breakout, we allow that the ejecta 
may no longer absorb the majority of the photoionizing 
radiation. We characterize this via the optical depth, 
through the ejecta at frequency z/, 

Tvif) = api{u,Xi) r, (7) 
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where api{v,Xi) is the photoionization cross-section and 
the sum is over all isotopes in X. The time an isotope can 
spend in the neutral state is set by the photoionization 
timescale, which we define at the ejecta outer radius (r = 
Veit), giving 


1 


( ^WD 


B,^(Tes) 

dv —--CTpi 

hv 


i^), (8) 


I^i/h 


where B,^{Tee) is the Planck function evaluated for the 
WD Teff, and Ixi is the ionization energy for isotope Xi. 

We consider only the two most significant elements for 
photoionization heating; H and He. We assume that after 
the breakout time n{Xi) is low enough that the ejecta is 
optically thin for Lyman limit photons {hv = 13.6 eV). 
In this limit we neglect any s econdary ionization ca used 
by a recombination emission ( Baker fc Menzel||1938 ) and 
consider only hydrogenic ionization processes (essentially 
assuming there are no neutral helium atoms). This gives 
two dynamical expressions for the hydrogenic number 
density of hydrogen and helium, respectively 


i{H) = - 


»(H) 

tion(H) 


-knen(H+)a(r,H+) ; 


(9) 


h(He''“) = —-I-nen(He^''')Q;(T, He'*”'') . (10) 


tio„(He+) 


As long as the WD is emitting, we found that the 
timescale over which these quantities relax is orders of 
magnitude smaller than the age of the ejecta, so we can 
consider an instantaneous steady state solution, where 
n{Xi) = Q, yielding 


n(W) — + Pci(T,X,) = nen(X+i)a(T, A+i) , 

\ ^ion / 

( 11 ) 

with Pci{T,Xi) the temperature-de pendent col lisional 
ionization rate which we adopt from Lang (1980). This 
equation yields the very small neutral fraction that we 
must know to calculate the heating rate from photoion¬ 
ization. 


3. TIME-DEPENDENT ENTROPY EVOLUTION 
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Figure 1. The calcu lated time depend ent Teff’s of different mass 
WDs as simulated by|Wolf et al.| (|2013|l, including two new cases 
(0.80 Sz 1.10 Mq) here, ihe age is the time since the end of mass 
loss driven by the thermonuclear instability. The rapid decline at 
later times for the three most massive WDs corresponds to the end 
of their SSS phases. 


with i?io„(Aj) being the average photoelectron energy. 
Under the op tically thin assumption Pion is given by 


Draine (2011) as 


Eion{Xi) — 


I dv ^^^^p-(jp-,{v,Xi){hv - 

Z^Iuth 

OO " ■ 

/ diy^^api{v,Xi) 

Z^Ia/h 


^ ^ (14) 

Using equations (1^ and (121 we hnd the equation for 
T(t) as 

dT 1 2T 

[rpi-A(r)]- —. (15) 


dt SkUe 


t 


Now that we have shown that photoionization balance 
drives the neutral fraction to a low value after the break¬ 
out of the Strbmgren sphere, we can consider the entropy 
evolution of the expanding ejecta, which will lead us to 
a temperature evolution. The photoionization heating 
timescale is much less than the age initially after break¬ 
out so that any information regarding the initial entropy 
is lost. This will simplify our calculation as we con¬ 
sider the state of the plasma under photoionization heat¬ 
ing and radiative cooling, where we use a cooling func¬ 
tion A(T) determined by comparative simulations run in 
Cloudy as detailed in Section The resulting rate of 
change of specific entropy s is then 

= Tpi - A(r) , (12) 

where the rate of photoionization heating, Tpi, is given 

by 


Finding a temperature evolution, T(t), requires solving 

equations (13) & (15) simultaneously. _ 

Before showing tne detailed solutions of equation (15), 
it is best to first exhibit a few limits. The first limit is 
to assume that the plasma is nearly fully ion ized, and 
that ionization balance is set by equation (11) but with 
little impact from collisional ionization. Tms yields a 
simple relation of nn/tion ~ UeUpU and a heating rate of 

<C t when the WD is hot, we 


Tn 


EionneiT-pa. Since tu 


found that ejecta of different initial temperatures would 
very rapidly (on a timescale much less than the age) reach 
a state of Tpi « A(r) and a temperature, Tiso nearly in¬ 
dependent of the electron density. The scale of these 
heati ng a nd cooling terms (which nearly cancel in equa¬ 
tion (15l) are so much larger than the 2T/t adiabatic 
expansion term that the ejecta evolves from one isother¬ 
mal (heating balances cooling) state to the next under a 
condition of thermal and photoionization balance. The 
resulting nearly isothermal evolution of this phase then 
arises from the condition of Tpi « A(T) for different WD 
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Teff’s. There is nearly no sensitivity to the actual elec¬ 
tron density when in this limit, so ejecta clumping will 
not matter, and all initial entropy information is l ost. 

However, the heating term, Tpi/rie, in equation (151 is 
oc implying that there will come a time, which we 
call tiso, beyond which the heating term is comparable 
to, or less than, the adiabatic expansion term, 2T/t, and 
we expect to see a temperature decline. If we neglect 
radiative cooling, this critical timescale, tiso, is when 
Tpi/Sfcrie = Eionnpa/3k = 2T/t. For the case of pure, 
ionized Hydrogen this gives 

1/2 , S -3/2 




250 d 


AMej 


10-5 

1/2 


Mq 

kT 


10 ^ km/s 


\10eV/ \1 eV/ V10-15 cm^ s-i / 

(16) 

We will show in the following section that our detailed 
thermal evolution calculations for these remnants can be 
simply described by the arguments above, with nearly all 
solutions exhibiting a nearly isothermal phase for t < tiso, 
followed by a temperature decline. 

4. CLOUDY SIMULATIONS 

We also simulated the nova eject a using the radiativ e 
transfer code Cloudy (version 13.03, Ferland et al.|2013). 



Fractional Mass Coordinate 

Figure 2. Temperature profiles of Cloudy simulations of nova 
ejecta with varying density profiles and compositions. Each repre¬ 
sents IQ-® Mq of ejected material with an outer velocity of 1000 
km s-^ being heated by radiation from a 1.00 Mq WD SSS. The 
abscissa measures the total enclosed mass at a given shell radius, 
with the outermost point being the position enclosing all the mass 


past by Schwarz 

(2002 

), Vanlandingham et al. 

(2005 

I, 

pchwarz et al. ( 

M7) 

Cloudy allows a user to spec- incorporate dynamic ettects other than the artihcially in 

he density structure ot a gas, its composition, and an serted adiabatic cooling. 


incident radiation field. It then computes the ionization 
state of the material, yielding line strengths, tempera¬ 
tures, and many other details, including realistic cooling 
curves. 

To make the models in Cloudy as similar to our semi- 
analytic approach as possible, we set the radiation field 
at a given time to be a blackbody spectrum with the 
Teff and radius set by the MESA models for that WD 
mass. The constant density gas reaches an outer radius 
Toutit) = Vejt and an inner radius rin(f) = O.Olrout(i) es¬ 
sentially reproducing our constant density sphere model. 
To account for temperature decline from adiabatic ex¬ 
pansion, we insert a “cooling” term 


Aoxp « 2 X 10 AMej _5 I^ej.3 ^4 tg ^ erg s 


cm-5, 


(17) 


that matches the second term in equation (|3), where 
AMej _5 = AMej/ (10-5 Mq), ■Uej,3 = i^ej/ (1000 km s ’' 


), 

T 4 = r/(10"‘ K), and = t/(10® s). The ejecta material 
was set to have a solar (Cloudy’s default abundance set) 
composition so that cooling by metals would be included. 
We then divide the WD evolution histories into twenty- 
one times, divided equally in log age, and use Cloudy 
to compute a temperature/ionization structure for the 
ejecta. These calculations assume that the cloud has 
reached both a photoionization and thermal equilibrium 
with the applied radiation field. We already established 
in Section that these are valid assumptions at early 
times while tion ^ t, so these Cloudy simulations pro¬ 
vide physically motivated snapshots at those early times 
of concern for us. Esse ntially, these simulations in their 
current form solve (15) assuming dT/dt = 0 and do not 


For these simulations, we found that the temperature 
was rather constant with radius, leading to a single tem¬ 
perature for each model that we defined as the mass- 
average over the outer 50% of the mass. These tempera¬ 
ture histories are what we compare to the semi-analytic 
models in Section For later use in our semi-analytic 
work (e.g. equation ( [T^ ) we also extracted a cooling 
curve by finding the volumetric cooling rate at each grid 
point, subtracting off the p orti on of that cooling due 
to expansion from equation 0 - and then dividing by 
riHU-e. The resulting cooling curves depend on the WD 
Teff (and hence WD mass) due to different degrees of ion¬ 
ization and line cooling. We thus generated characteristic 
Cloudy cooling curves for each WD mass by choosing a 
time in the Cloudy calculations about halfway through 
the observed temperature plateau in log age. At times 
beyond these, the heating and cooling timescales become 
comparable to the age, and we can no longer trust the 
equilibrium Cloudy calculation, but rather should solve 
the fully time-dependent problem. We make progress on 
that aspect in the Section 

Since Cloudy is inherently a multi-zone code, we can 
use it to test our assumption that the density profile 
and composition has little influence on the temperature 
structure of the ejecta. We ran four models with our 
fiducial ejecta parameters of AMej = 10“® Mq and 
ifej = 1000 kms-^ with the radiation field given by the 
1.00 Mq WD model. These four models used two dif¬ 
ferent density distributions, constant and inverse square, 
as well as two sets of elemental abundance ratios, solar 
and nova. The nova abundances are a built-in option 
in Cloudy and are intended to model the abundanes ob- 
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served in the ejecta of nova V1500 Cygni (Ferland & 


Shields 1978). Since the adiabatic cooling is necessar- 
ily radius-dependent for a non-cons t ant density profile, 
the simple extra cooling of equation 0 could not prop¬ 
erly model expansive cooling. To be consistent, we in¬ 
cluded no extra cooling from expansion in this compar- 
ision study, resulting in hotter temperatures and longer 
“plateau” periods. The temperature prohles of all four 
cases at a time approximately 200 days after ejection are 
shown in Figure Importantly, the temperature range 
changes no more than by order unity in the constant den¬ 
sity models and is nearly exactly isothermal in the inverse 
square cases. Also, the temperatures for these models do 
not change greatly as the density profile and/or com¬ 
position are changed, allowing us to use our simple so¬ 
lar composition and constant density assumptions in our 
semi-analytic model. One important difference between 
the different density profiles is that the opacity at the 
lyman-limit diverges as we let the inner radius go to zero 
in the inverse square case. While a vanishing central ra¬ 
dius is certainly unphysical, we must be aware that the 
time to become optically thick is sensitive to the choice 
for an inner radius. In this comparison we continued to 
use Tin = O.OlUejt, which causes the initial optically thick 
phase to last significantly longer than the similar con¬ 
stant density models, sometimes approaching 100 days. 

5. RESULTS 

The results of the one zone semi-analytic model and 
Cloudy simulations for WD mass 0.60, 0.80, 1.00, and 
1.10 Mq are shown in Figure and Table from the 
time-dependent WD models described earlier. All tem¬ 
peratures are in the 10“^ K region, consistent with that in¬ 
ferred from radio observations and nebular spectroscopy, 
and show a long plateau of nearly constant temperature 
at the value expected for photoionization and thermal 
equilibrium. Results for the two higher-mass WD mod¬ 
els are not shown for reasons described below. 

The four WD evolutions shown exhibit a period of 
isothermal expansion lasting around 100-200 days that 
we will refer to as the plateau. Numerically, we de¬ 
fined this as the duration over which the the tempera¬ 
ture stayed within 20% of the plateau temperature Tiso- 
These durations are i^ accord with the timescale pre¬ 
dicted by Equation (16), and show that, even with the 
WD still providing a nrgh level of photoionization, the 
adiabatic expansion term eventually plays an important 
role. 

The two methods - one-zone model and Cloudy simu¬ 
lations - agree well for all four cases. The two higher 
WD mass simulations (1.20 Mq and 1.30 Mq) showed 
much shorter plateau phases (up to 80 days). However, 
this was due to the sharp decrease in Teff at early times 
relative to the 4 lower mass WDs, see Figure [H The 
plateau phases were so short in these cases that the 
ejecta remained optically thick at the photoionization 
edge throughout and in our model these solutions would 
not have been self-consistent. For this reason we consider 
only the WDs in the mass range 0.60 — 1.10 Mq. Note 
that for lower ejecta masses and higher ejecta velocities 
expected of novae on higher-mass WDs, it is possible 
that the plateau phase would occur after the ejecta has 
become optically thin. In this work, though, we hold 
the ejecta parameters constant and thus do not consider 


Table 2 

Isothermal Evolution of Nova Ejecta** 


Myjy) 

(Mq) 

q-' b 
iso 

(lO'^ K) 

f- ^ 

‘'ISO 

(d) 

hhin'^ 

(d) 

0.60 

2.2 

229 

36 

0.80 

3.2 

155 

33 

1.00 

4.7 

152 

47 

1.10 

4.5 

217 

74 


** AALej = 10 ® Mq; Vgj = 10® kms ® 

** Time-averaged temperature in isothermal phase 
** Time after temperature drops 20% from early times. 
Duration of optically thick (to lyman-limit photons) phase 


the resulting inconsistent results from higher-mass WDs. 
Most of these results can be motivated and explained by 
our analytic work that derived equation (16), yielding the 
dependence of tiso on the ejecta parameters and the aver¬ 
age photoelectron energy. Such an understanding will be 
key to using these observed quantities to constrain either 
the WD mass or the ejecta properties. 

To exhibit this physics, we explored the sensitivity of 
the plateau phase properties (Tiso and tiso) by varying 
Uej, AMgj and Aion in simulations of a pure H ejecta pho- 
toionized by a 1.00 Mq WD. We found that variations 
in ejecta properties (ugj and AMgj) affect the duration of 
the plateau but do not alter its temperature. This is re¬ 
vealed in Figure 15 where the ejecta mass varies between 
10“® Mq < AA/^ < 10“^ Mq with velocities in range 
10^ kms“^ < I'ej < 5 X 10^ kms“^. For calculating Fig¬ 
ure 15 the Aion was held constant in time for simplicity. 

This highlights that the dynamic behavior of the ejecta 
is described by these two quantities, Ugj and AMgj, that 
we can parametrize into a factor 




AMg. 

10“^ Mq VlO® kms“® 


( ) 
V 10 ^ kins“^ / 


-3 


(18) 


Using Equation (16), ( can be inferred from the observed 
parameters tiso ana Tjso as 




potentially yielding insights into the ejecta properties 
from observations. The average photo-electron energy 
Eion depends on the WD’s Tgff and is the main determi¬ 
nant of Tiso, as is evident in Figure]^ Hence, the plateau 
temperature is most sensitive to the WD Tgff (and hence 
WD mass). 

The plateau phase space behavior of is shown in Fig- 
ure|^with dashed lines of constant ((. For nearly all cases, 
we see that Tiso is highly independent of the ejecta infor¬ 
mation, Ugj and AMgj. This implies that changes made to 
the dynamic properties of the ejecta should be relatively 
insignificant for the temperature, making it a potentially 
clean diagnostic for WD mass. This also implies that 
changes to the ejecta profile (e.g. non-spherically sym¬ 
metric ejections) or homogeneity would unlikely impact 
T 

ISO • 

The thermal evolution of the ejecta at late times is sim¬ 
ple if the SSS phase has ended, as the ejecta undergoes 


Ei 

0 


eV J V] 


Q 


r), (19) 
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Age (days) 


Age (days) 


Figure 3. Temperature evolution of the nova ejecta for four WD masses using our semi-analytic model (blue) with cooling curves extracted 
from Cloudy and simulations run in Cloudy (dashed gold). AMgj = 10“^ Mq and Uej = 1000 kms“^ for all models, and the time evolution 
of the incident radiation is set by the underlying WD mass via the histories shown in Figure^ The dotted line indicates that the ejecta is 
optically thick at the Lyman limit edge at early times. 


an adiabatic temperature decline, and T oc there¬ 
after. The recombination time is far longer than the age 
when t > tiso, so that this plasma would maintain itself 
in a highly ionized state even though photoionization has 
ended. Hence, it would still be a thermal bremsstrahlung 
radio emitter during the adiabatic decline. If the SSS 
remains on at t > tiso, we can still approximate the ther¬ 
mal evolution. In this limit, photoionization heating is 
still operative, so the temperature decline would not be 
as steep as the adiabatic relation. Rather, for the cases 
we considered here, the photoionization timescale is still 
much less than the age, tion ^ t, so that the heating re¬ 
mains operative, and the required small neutral fraction 
is easily maintained. If we assume negligible radiative 
cooling, rather safe as T < K, then equation (151 
simplifies to an integrable form. If the ejecta tempera¬ 
ture was To at t = to, then its value at a time t > tjso is 
simply 


T 


t2 




'ISO 

^0 


t- 

^ISQ 


( 20 ) 


exhibiting the less steep than adiabatic decline that is 
evident in nearly all our results. Differentiating be¬ 
tween a decline from SSS shutoff versus onset of dom¬ 
inance of adiabatic expansion with an active photoion¬ 
ization source could be attempted by comparing T{t) in¬ 
ferred from radio observations or nebular spectroscopy to 


T oc versus the more complicated relation of equa¬ 
tion (20). Alternatively, a shut-off of the X-ray source 
could be confirmed by direct measurements like those of 


Schwarz et al. 

(2011 

) and IHenze et al. ( 

2011 2014b 


the WD rnodels used in Figure are still increasing 
in effective temperature when adiabatic expansion t akes 
over, so their decline is even less steep then equation (201 
would indicate. 


6. CONCLUSIONS 

We have shown that the photoionizing emission from 
the hot white dwarf following a classical nova is ade¬ 
quate to explain the inferred ejecta temperatures > 10"^ 
K that persist for nearly a year after the optical out¬ 
burst. For the typical WD mass, ~ O.SMq, of a CNe 
this nearly isothermal phase lasts for about one year (see 
equation (16)) and has a temperature of T « (2 — 3) x lO'^ 
K, in good agreement with many observed novae. Even 
when the SSS phase lasts for decades, the balance be¬ 
tween heating and cooling eventually ends, and the ejecta 
transitions into a phase of temperature dec line roughly 
approximated analytically by equation (20). If the SSS 


phase ends early, then the ejecta temperature would 
evolve adiabatically, as T 

Though our preliminary calculations allow us to see 
the basics of how photoionization heating is balanced 
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Figure 4. Ejecta temperature evolution for variations in ejecta 
velocity (top panel) and ejecta mass (bottom panel). Simulations 
carried out assuming pure H ejecta with constant Eion and neg¬ 
ligible contribution from collisional ionizations (P^i = 0)- Unless 
otherwise indicated AMej = 10“® Mq, Uej = 10^ km s“^ and 
Eion = 22 eV. 


by cooling at early times, and how the heating is weak 
compared to adiabatic expansive effects at late times, 
work remains to realize the quantitative accuracy ade¬ 
quate to infer CNe ejecta properties or WD masses. We 
certainly identified a degeneracy in how the ejecta infor¬ 
mation (AMej and Uej) affect the plateau phase duration, 
which should allow observers to make predictions about 
the properties of a given nova, complementing existing 
rad io techniques dSeaquist fc Pa limaka||1977[ |Hjellming 
d.|1979 Heywo^ et al.|20M). h'or example, equation 


et al.| 


1 19) makes possible the inference of an ejecta mass range 
based on the length of the isothermal temperature evo¬ 
lution inferred in the radio combined with the measured 
velocity of the ejecta and temperature. 

However, specific events would be best modeled by a 
unique Cloudy simulation that reflects the known proper¬ 
ties of the event as inferred from the optical, as well as the 
X-ray supersoft phase. That would also allow for inclu¬ 
sion of unusual abundances c onstrained from modeling 
the nebular spectrosc opy a s in Schwarz (12002 b. Vanland- 


ingham et al. (2005), and Schwarz et ^12007), which 




01 

d 

01 
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Figure 5. Ejecta temperature evolution for varying average en¬ 
ergies of photoelectrons from 4 eV < Eion < 22 eV. Simulations 
carried out assuming pure H ejecta with constant -Eion and neg¬ 
ligible contribution from collisional ionizations (Pci = 0)- Unless 
otherwise indicated A.Mej = 10 ~^ -^©j = 10 ^ km s~^. 
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Figure 6. Characteristic plot of the plateau phase with dotted 
lines of constant ejecta parameter = 10 ^"^ x — 3 -^ spaced evenly 

in log space. A line of a particular color connects models of identical 
WD mass. 


would in turn impact the mapping of the observed ejecta 
temperature to the WD mass. The ideal calcuation could 
use a more realistic ejecta dens ity profile, like a Hubble 
flow (Seaquist & Palimaka|1977 ) with an appropriate fill¬ 
ing factor, though our efforts irnply that the ejecta tem¬ 
perature is relatively independent of the specific density 
structure for the first year or so. We hope that the ex- 
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panded efforts with the Karl G. Jansky Very Large Array 
to monitor nearby galactic novae will enable a new era of 
quantitative ana lysis of the late time radio observations 
(|Roy et al.||2012|). 
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